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ABSTRACT 

The behavior of correlation functions is studied in a class of matrix models characterized by 
a measure exp(— S) containing a potential term and an external source term: S = N tr( V(M) — 
MA). In the large N limit, the short- distance behavior is found to be identical to the one 
obtained in previously studied matrix models, thus extending the universality of the level- 
spacing distribution. The calculation of correlation functions involves (finite N) determinant 
formulae, reducing the problem to the large N asymptotic analysis of a single kernel K. This 
is performed by an appropriate matrix integral formulation of K. Multi-matrix generalizations 
of these results are discussed. 
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1. Introduction. 



More than four decades ago, Wigner [1] suggested to study the distribution of energy levels 
of complex systems using random matrices. In this approach, one would like to characterize 
the structure of the energy levels, i.e. the eigenvalues of the Hamiltonian, the latter being 
considered as a large matrix with random entries. It is now known that many statistical 
properties of spectra of true physical systems are indeed well described by those of random 
matrices (cf [12] for a review): it is therefore important to understand how much these spectral 
properties depend on the particular matrix ensemble chosen, i.e. determine universality classes 
of matrix ensembles. 

For technical reasons, we shall consider here ensembles of hermitian matrices, which cor- 
respond to systems without time-reversal invariance. The main quantities of interest are the 
probability distributions p n of the eigenvalues: if M is a random hermitian N x N matrix, 
we define p n (Xi, . . . , X n ) to be the density of probability that M has (Ai, . . . , A n ) among its 
N eigenvalues, with the normalization convention that: J* Y\H=\ dXi p n (Xi, . . . , A n ) = 1. To 
connect with correlations functions of the model, we define (following [2]): 



coinciding eigenvalues). This means that correlation functions of any U (iV)-invariant quantities 
(i.e. functions of the eigenvalues only) can be computed using the p n . 

We shall now study a class of models in which one can express the functions p n in terms 
of a single kernel K(X, p); if we again assume the X{ all distinct, the corresponding relation for 
FL will be 




(1.1) 



Then one has R n (Xi, . . . , X n ) 



p n (X\, . . . , A n ) for distinct Aj (5 functions appear for 



{N-n)l 



Rn(Xi, Ai, . . . , A n ) = det(K(Xi, Xj))ij=i... n . 




These formulae are exact at finite N. 



For example if we also define 




(1.3) 



l 



the connected correlation functions, this implies that 

R i 2 c \\, f i) = -K(\, f i)K(fi,\). (1-4) 

In these models, the study of the distribution of eigenvalues reduces to the analysis of this 
kernel; in particular the large N limit of K should allow to compute all correlation functions 
in this limit and find the different "universal" behaviors that can arise: 

* The long distance behavior. As it is known that the kernel fluctuates wildly on intervals 
of size ~ 1/N, one must first average the kernel to suppress the oscillations on this scale and 
obtain a sensible long distance behavior. 

* The short distance behavior. This is the region X—pL ~ 1/iV, in which the fast oscillations 
mentioned above are relevant. 

The long distance behavior can be studied by various standard large N techniques [7,11,14], 
so we shall concentrate here on the short distance behavior. Usually one characterizes this 
behavior by introducing the level-spacing distribution P(s), s = N(\ — jj) (P(s) a priori also 
depends on A or /i). In the large N limit, P(s) can be simply related to the asymptotic form 
K of the kernel on [A, //]: 

(1.5) 

where det is the Fredholm determinant. Thus, in the short distance region, the universality of 
K implies the universality of the level spacing. 

In section 2 we first consider the case of a matrix model, where the measure consists of a 
simple potential term tr V(M). Expressions of the kernel K in terms of orthogonal polynomials 
have been known for a long time [2] . One can then proceed to derive a short distance universal 
behavior of the kernel [11,15]: 

si n t 

X(A,aO~-— , x~N(\-ii). (1.6) 

Ob 

Here we shall rewrite K, and rederive its large N asymptotic form, using a new method which 
does not make use of orthogonal polynomials, keeping in mind that we are ultimately interested 
in the more difficult case of the model of matrices coupled to an external field. 



P(s) 



ds' 



det 



1 K 

N 



2 



The first hint that the latter model could possess the same short distance universal behavior 
appeared in a series of papers [17] by Brezin and Hikami, who showed that formula (1.2) 
can be generalized to the gaussian ensemble with an external field, that is for the measure 
exp (-f trM 2 + Nti MA) d N M . This measure can be interpreted by saying that M is the 
sum of a fixed Hamiltonian A and a random gaussian part M — A: the constant part A breaks 
the U(N) invariance of the model so neither orthogonal polynomials (like in the standard one- 
matrix model) nor biorthogonal polynomials (like in the two-matrix model) are of any use. 
Still, one can define a kernel in this model, which has the same short distance behavior (1.6). 

Then in [18] the more general case of a measure of the type exp(iVtr(— V(M) + MA)) 
was investigated. The motivation of this measure was to study a random Hamiltonian which 
contains a not necessarily gaussian potential, and an external source term which breaks the 
£/(iV)-invariance: could the latter symmetry be somehow related to the universality of the 
level spacing distribution ? 

To answer this question, we shall again define in section 3 a kernel K such that equation 
(1.2) holds. Then, just as in section 2, we shall rewrite K as a matrix integral, allowing to 
take the N — > oo limit, and find the short distance universal behavior. 

The methods used here are very general, and in particular, they can be successfully applied 
to multi-matrix models; to show this we study in section 4 a model of a chain of matrices, with 
or without an external field at the end of the chain. We give without proof, and in analogy 
with the one-matrix model, expressions for the kernel K, which exhibit the same short distance 
universality 

Finally, appendix 1 describes in detail the analytic structure of the functions involved in 
the large N limit, and appendix 2 shows the connection between the formalism used here and 
large iV character formulae. 
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2. The U (iV)-invariant case. 



2.1. Definition of the model and of the kernel. 

Let us consider an ensemble of random hermitian N x N matrices with the measure 

Z _1 exp(-iV tr V(M))d N2 M (2.1) 

where V is a polynomial and Z the partition function. An important remark is that this is 
not the most general U (AQ -invariant measure (one could have products of traces of functions 
of M). 

A classical result [2] expresses the distribution law p n of n eigenvalues (1 < n < N) of M 
in terms of the kernel 

N-l 

K(\,n) = J2 F k( X ) F M- ( 2 - 2 ) 

fc=0 

Here F{ is the orthonormal function associated to the usual orthogonal polynomial -Pi (A) = 
A* + ---: 

Fi(\) = h-^ Pi {X)e-^) 

(2.3) 



dXe-^WPiWPjiX) = hid. 



ij- 



(see [15] for a review of orthogonal polynomials in matrix models). Let us briefly rederive this 
result. As the measure (2.1) only depends on the eigenvalues of M, the integration over the 
angular variables is trivial and one finds: 



(Ax, A 2 , . . . , X N ) = Z-^AOe""^! V(X < (2.4) 



PN 



The Van der Monde determinant A(Aj) = det(Ai- ? )i<j<jv ) o<j<iV-i can De rewritten in terms of 
the orthogonal polynomials: 

Piv(Ai, A 2 , . . . , Ajv) = Z" 1 det(Pk(\i))i<i<N,o<k<N~i det(P k (\j))i<j<N,o<k<N-ie~ N ^i=i y(Al) . 

(2.5) 

One can now easily compute Z = N\ n^Q 1 ^* by integrating over all A,. Combining the two 
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determinants, we finally obtain: 

Piv(Ai, A 2 , . . . , Ajv) = det(K(Xi, \j))ij=i...N- (2.6) 
The kernel K has the following properties: 

f (2.7) 
[K*K\(\,v)= dnK(X,n)K(n,u) = K(X,u) 

i.e. it is the orthogonal projector on the subspace spanned by the F^, < k < N — 1 . Using 
the property K * K = K and noting that 

Pn(Al, A2, • • • , A n ) = J dAn+lPn+l(Al, A2, . . . , A n+ i), (2.8) 

one can then show inductively that 

Pn(Ai, A 2 , . . • , An) = V m ! ' det(K(A i? A,));.,- ,...„. (2.9) 
for any n < N . This is equivalent to formula (1.2). 
2.2. Matrix integral formulation of the kernel. 
We shall now rewrite K(\,/i) as a matrix integral: 

K(\,n) ^Z~ 1 e-^^ +v ^ J S- N ~^M det(A - M)det(/i- M)exp(-NtrV(M)). (2.10) 
It is easy to check this formula by going over to eigenvalue variables: 
J rf (iV " 1)2 Mdet(A - M) det(/i - M) exp(—NtrV(M)) 

f N-l 

= / ]J(dX l (X-X l )(fi-X l ))A 2 (X l ) l < i < N ^e~ N ^ lv ^ 

J i=\ 

.N-l _ (2.11) 

= / n^ A ( A *)i<^A Ar .AA(A0i<,<^^e^^i 1 ^) 

i=l 

N ~ l r r 

= E (-in-l)' P*mWPo>(N)M II J dzP a{i) (z)P af{l) (z)e- NV ^ 



This is non-zero when a = a' and we find as expected: 



N-l 

K(X, = e-T(VW+VM) h- l P k {\)P k {n). (2.12) 

fc=0 

Note that equation (2.10) does not involve orthogonal polynomials, and that is why we shall 
be able to generalize it to the non £/(iV)-invariant model of matrices coupled to an external 
field, for which the orthogonal polynomials formalism is not available. 

2.3. Large N asymptotics of the kernel. 

Formula (2.10) allows to compute asymptotics of K(X,fx) as iV — > oo. It relates K(X,/i) 
to the partition function Z(X,fj) of a matrix model with the measure exp(trlog(A — M) + 
trlog(/i — M) — NtrV(M))d( N -V 2 M: 

Z(\,n) = y"rf (iV " 1)2 Mdet(A-M)det(/i-M)exp(-iVtr\/(M)). (2.13) 

Rather than directly applying the saddle point method to this expression, it is easier to write 
differential equations for Z. Indeed one has: 

^-logZ(A,^) = (JV-l)G A>/1 (A) 

* (2.14) 
— logZ(A,/x) = (iV-l)G Ai » 

where is the resolvent of this model: 

G x Jz) = — — -(tr — l — ) (2.15) 
; N- 1 \ z-M I y ' 

which depends on A and /x through the trlog(A — M) and trlog(/i — M) terms in the action. 
Note that the factor iV — 1 in front of G\ yfI in (2.14) forces us to compute G up to l/N 
corrections. 
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Let us now find a saddle point for the eigenvalues. In the large N limit, we shall suppose 
that they fill a single interval [a, /?]; then Gx^(z) becomes an analytic function of z with a 
single cut on [a, j3\: 

G x ^z ± iO) = fixM ± ™P\A Z ) ( 2 - 16 ) 

where p\^(z) is the density of eigenvalues at z. As a definition G(z) = ^(G(z + iO) + G(z — iO)) 
for any function G. The saddle point equation can now be written: 



2(N - l)GxM ~ NV '( Z ) + ~T + — = Vz e [a, 6]. 

z — A z — ji 



(2.17) 



At leading order in JV, this equation is just the usual saddle point equation 2G(z) = V'(z) for 
the resolvent G of the original matrix model we have started from. Therefore we can write: 



Gx,» = G+-(C G + Cx + C fl ), 



(2.18) 



where we have introduced three 1/N corrections to the leading behavior of Gx,^ corresponding 
to the three corrective terms in the saddle point equation (2.17). Following [10,15] we deduce 
these corrections from their analytic properties. Both G(z) and Gx^{z) behave as \jz as 
z — > oo, so Co, Cx and C M are 0(1/ z 2 ). Furthermore they satisfy 



f <p G (z) = G(z) 



<P»( Z ) 



2(X-z) 
1 

2(A* - z) 



(2.19) 



and are regular for all z except near a (resp. (3) where they should behave as 1/y/z — a (resp. 
1 / y/z — (3) . This determines them entirely: 



C G {z) = G(z) - 
Cx(z) 1 



y/(z-a)(z-(3) 
1 ( ^ y/{z -a)(z-(3)- y/(A -g)(\-/3) 



C,(z) 



2 v/^-a)^-/?) \ * 

I f i y/(z- a)(z -(3)- y/(/i -a){n-l3) 



(2.20) 



2 V(* -«)(*-£) 



z — n 
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The right hand side of (2.14) can now be written: 



(2.21) 

and a similar equation for (N — l)G\^(n). For A G [ct,/3], an ambiguity in (2.21) must be 
resolved: G\ yfI , like G, has a cut on [a, /3]; so we must choose A slightly above or below the 
real axis to determine the right hand side of (2.21): 

(N - l)G Arfl (A ± iO) = ^V(X) ± Ninp(X) - ~ log ^(A -«)(/?- A) 

1 1 ( 1 V(pt-a)(pt-P) N (2 ' 22) 

2 A - /i y ±2a/(A - «)(/? - A) 

(p(X) = pi(A)). This ambiguity, which appears only at iV = oo, means that there are several 
saddle points which we must all take into account. The same problem appears when p, gets 
close to the cut, so there is a total of 4 saddle points (e, e') (e, e' = ±) corresponding to the 
locations of A and p with respect to the cut [a,/?]. 

Finally we can write differential equations for K(X,p): 
— XogK^ = dVmp(A) - -— log a/(A - a)((3 - A) - -— - 

9 . „ , s lrf, „ — r 1 1 



a) sin<£>(;z) = ^/(z — a)(/3 — i). Noting then that 





«)(/?- 


**)\ 


ev /(A- 


«)(/?- 




ev /(A- 


-a)(/3- 






-a)(/3- 








(2.23) 


ce) cos<pl 


'2) and 





d logGin Ay(A)-6V(/i) A _ c 1 ev /(A - a)((3 - A) + e'^ -«)(/?- //) 



dA ° V 2 ; a/(A - a) (j3 - A) A — // 

(2.24) 

equations (2.23) can be integrated: 



sin I 2 J 1 if if 

K {t,e')(\v) = C (e,e') T / ■ 7TY=- =T^ eX P diV7r / P(z)dz + e'lNlT / 

A — /J ysm </?(A) sm </?(/j) \ J J 

\ Ao Ao 

(2.25) 

The integration constants C( ee /) satisfy C( ee /) = C( £ / e ) (interchange of A and p) and C(_ e _ e /) = 
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— ci ee i\ (complex conjugation). c {±,^) are independent of the choice of Ao and are fixed by 
imposing the normalization condition K(X, A) = Np(X): we find (cf next section) that cr± ^ = 
l/{2m). C(±,±) are undetermined; if we assume that one can find Ao such that C(±^±\ = ±l/2n 
(for the case of an even potential we would have Ao = 0), then we are left with only one 
unknown parameter Ao- 

We can finally sum the four function Ku e i\] we obtain: 
1 1 1 



K A — fi sin (/9(A) sin(^(/i) 

sin( m±M ]siIi {NlT j P ( Z )dz 




(2.26) 



+ sin ^(A) _ v(fJ.) \ cog ( Njr J p ^ dz + Nn J p ^ dz 



a formula for K that is equivalent to the one that was found in [11] by using an ansatz on the 
form of orthogonal polynomials (see also [15]). 

2.4. Short distance universal behavior of the kernel. 

Let us now inspect the region where A — /i ~ 1/iV, a < A < (3. It is clear from (2.25)-(2.26) 
that the dominant contributions come here from the saddle points (±, =f) (i.e. A and /i on 
opposite sides of the cut). Actually, this can already be seen in the equations (2.23), which 
acquire a particularly simple form in this limit: 

, . 1 ^ 



Here we do not need these simplified differential equations, since we can directly take the limit 
in (2.26), which yields 

K(A, / ,) = Sin(iV7r ^-^ (A)) . (2.28) 
This is the well-known short distance universal behavior (1.6) of the kernel. 
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3. Generalization to the case of an external field. 



It was proven in [18] that in the case of a general measure with an external field, (1.2) 
still holds; we shall now review this result and write the kernel K in an appropriate way for 
asymptotic analysis. 

3.1. Definition of the model and of the kernel. 
Let us consider the measure: 



where V is an arbitrary polynomial, and A = diag(ai, . . . , ajv) can be assumed diagonal. 
Particular matrix models of this type appear in several papers [6,8]. 



Z' 1 exp (-N trV(M) + Ntr MA) d N M 



(3-1) 



One diagonalizes M: MM — VtAVfi where A = diag(Ai, . . . , Ajy), the integral over VL is the 
usual Itzykson-Zuber integral [4] on the unitary group and we find: 



Pn(M, A2, • • • , Ajv) = Z 1 A(X i ) 



det(exp(Aq,az)) -NY. N 1 v(\<) 



(3.2) 



A(o z ) 



Z can now be computed: 




(3.3) 



Inserting (3.3) into (3.2) yields 



pjv(Ai, A2, . . . , Ajv) 
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putting together the three determinants (and the exp — NV(X) factors) we finally obtain: 

Piv(Ai, A 2 , . . . , Ajv) = det(K(\i, Xj)) i;j= i__. N (3.5) 

where 

N-i N 

K(X, = e-T (VW+VM) J- a kl X k e Na ^. (3.6) 

k=0 1=1 

The kernel K satisfies the property: 

' ' ' dnn k 'e Na ^~ NV ^ 



[K*K](X,p)=e-T(V( x )+V( P )) 
= K(X,p). 



L.I 1 II \-J 



a kT e Na i' p 



(3.7) 



Thus, one can follow the same line of reasoning as in the (AQ -invariant case to obtain the 
determinant formulae 



(TV - n)\ 

Pn(Xi, A 2 , . . . , A n ) = — — — det(K(Xi, \j))i,j=i... n (3. 



for any n < N. 

If we introduce the polynomials Qf. 



N-l 

k 



QiW = £ aklX 

k=0 

then 

N 

K(X, /i) = e-T (VW+VM) J- Qi(X)e Nai ». 

l=i 

The polynomials Q\ are of degree N — 1, and satisfy the orthogonality relations: 

J rfAQKA)e iV( - y(A)+a; ' A) = 6 W , (3.9) 

This proves that if is a non-orthogonal projector on the space spanned by the 
Ql(X) exp(-W(A)/2), 1 < I < N, which is also the space spanned by the F k , < k < N - 1. 



11 



3.2. Matrix integral formulation of the kernel. 



We shall now guess matrix integral formulae for the polynomials Q[\ 

Q,(A) =c t J rf (iV " 1)2 M det(A - M) exp(Nti(-V(M) + MA®)). (3.10) 

Here A® stands for the diagonal (N — 1) x (N — 1) matrix obtained from A by removing the 
eigenvalue a\. The q are normalization constants. The right hand side of (3.10) is obviously 
a polynomial of degree N — 1. As the polynomials Qi are entirely characterized by property 
(3.9), we compute 

J d\iQi(\i)e N ^~ v ^ +aXl ^ 

i , N (3.11) 

= CI ( N - /n^ A ^)i<^^ E ^° ( " nAi)+aiA!KW(Ai)+ ^ 

where we have introduced the eigenvalues \, % ^ I, of the matrix M, and A(A^) = A(a//)//^/ 
is the Van der Monde determinant of the eigenvalues of A^. Eq. (3.11) looks like a N x N 
matrix integral; if we define A® ,a to be the diagonal N x N matrix obtained from A by 
replacing a\ with a, we have: 

d\ l Q l (\ l )e N (- v M +aX * = ql ^j^ y ^ 2 Mexp(iVtr(-K(M) + Mi^ ,a )). (3.12) 

If we now set a — a^, I' ^ I, the Van der Monde determinant A( y A^ ,a '') of the eigenvalues 
of A®*' becomes zero. If a = a/, A^)' a! = A and the matrix integral is just the partition 
function Z. For (3.9) to hold we need c\ to be: 

ci = NZ- 1 - i -. (3.13) 



(3.14) 



We can finally express the kernel as: 

r N-l 

i<i<N,x N =x det(exp(A r Aja^))i<j ) Kiv ) Aiv^- 

Note that itself, in contrast with the polynomials Qi, is more naturally expressed as an 
integral over eigenvalues than as a matrix integral. 
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3.3. Saddle point equation and analytic structure 



Before going on with the study of the kernel, we need to understand the analytic structure 
of the various functions that we shall now introduce. To do so, we first write saddle point 
equations for the standard partition function: 

r N 

Z~ / JJrfA i A(Ai)det(exp(A^A i a;))e- iV ^=i 1/ ( Al ). (3.15) 
"> i=i 

We shall suppose that the density of the N eigenvalues of A has a smooth limit as N — > oo. 
Then the eigenvalues of M also have a smooth large N limit, characterized by a saddle point 
distribution: we shall assume that the eigenvalues fill a single interval [a, (3} (it will be argued 
later that this hypothesis is only technical and does not change the short distance universal 
behavior), with a density p(A) = pi(A). 

Usually, at this stage, one replaces the determinant det(exp(iVAja/)) with exp(NJ2^i a i) 
using the symmetry of exchange of the eigenvalues; here we shall not do so, because this would 
prevent us from writing down a saddle point equation. Instead, we introduce 2 functions G(z) 
and a(z) which have a cut on such that: 

w 4^ logA(Aj) 

x Q ( 3 - 16 ) 
= — ^-logdet(exp(iVAia;)). 



G is of course the resolvent: 



a cannot be defined by such a simple formula; we refer to appendix 1 for a rigorous definition 
of a. Here let us note that Itzykson-Zuber formula 



det(exp(iVAia/)) 



J dQe Nt ^ m ' A \ (3.18) 



A(Ai) 

QeU(N) 

implies that d/d\i log(det(exp(A r Aja/))/A(Aj)) is a regular function of \ (i.e. there is no pole 
when Aj ~ Xj); so we write it under the form /(Aj), where / can be extended into an analytic 
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function on the whole complex plane. Thus a(z) and G(z) are related by: 

a(z) = G(z) + f(z). (3.19) 
From (3.19) we deduce that a has the same cut as G on [a, f3], that is 

a(z ± iO) = 4{z) ± iTrp(^) G [a, /3] (3.20) 

(3.20) can be taken as a definition of a on [a, 0\ since ^(z) is given by (3.16) (see appendix 1 for 
a definition of a(z) on the whole complex plane and a more detailed analysis of its properties). 

The saddle point equation for (3.15) is 

G{z) + i{z) = V'{z) Wz e [a, /3] (3.21) 

Using what we know of the analytic structure of G and a (eq. (3.20)), we can now extend 

(3.21) to the whole complex plane: 

G(z) + a*(z) = V'(z). (3.22) 

a* denotes the function connected to a by the cut [a, f3]. In other words a is a multi- valued 
function of z, and a(z) and a*(z) are the values on opposite sides of the cut since (for 

generic A and V) the cut is a square root-type cut which connects two sheets (figure 1). 




Fig. 1: Analytic structure of a(z) and G(z). In the physical sheet (below), G(z) 
has a single cut, while in the other sheet (above) it can have more cuts. It is the 
opposite for a(z), because of eq. (3.22) 
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3.4. Short distance asymptotics of the kernel. 



Let us analyze the more complicated saddle point equation for the integral 



N-l 



Z ( X ^) = / II ^ e "^ E ^ 1V(Al)A ( A Oi< l <^A Ar ^Adet(exp(ATA,aO)i< t ,Kiv,A J v^ (3.23) 



which is related to the kernel K by eq. (3.14). As in the £/(iV)-invariant case, Z(X,fi) can be 
considered as the partition function of a model in which the action is the action of (3.15) (of 
order N 2 ) plus additional terms dependent on A and /i (of order N). 

We shall again write differential equations for Z; in a very similar way to the ?7(iV)-invariant 
case, we find here: 



with a modified saddle point distribution of the A^ due to the additional terms in the action. 
Of course, the leading behaviors of G\^ and ax^, when N — > oo are simply G and a, since 
the corrective terms are negligible in the leading approximation. 

If we wanted to solve the differential equations for all values of A and fx, we should now 
calculate the 1/N corrections to the leading behavior. However, as we are only interested in the 
short distance behavior of K, we can restrict ourselves to the region A — fi ~ 1/N: variations of 
A, n around the diagonal A = ji (where K is known - K(\, A) = iVp(A)) are then of order 1/N, 
i.e. we only need the leading behavior of d/d\\ogK. This means that the next corrections of 
G\^(z) and a\^(z) are actually irrelevant in this region, except for possible poles of the type 
l/(z — A) or l/(z — //), which would be again of order N. 

Of course G\^ and do not have any poles on the "physical sheet" ; but we have learnt 
from the U (iV)-invariant case that different choices of sheets (or of values on the cut joining 
these sheets, which amounts to the same) correspond to different saddle points, and that is 
how we were led to taking into account saddle points (±, =f) in which poles at z = A and z = \i 
do appear (cf figures 2 and 3 in next section). 



i=l 




(3.24) 



G\^ and are defined in the same way as G and a (see appendix 1 for a definition of a), but 
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We now recall that a(z) = G(z) + f(z), where / is a an analytic function (regular on the 
cut [a, /?]). In the same way a\^(z) = G\^(z) + fx,u{z) and it is now clear that poles can 
only come from the Van der Monde part of a(z) and not from the regular part f(z). More 
explicitly, one can write down a saddle point equation for and which looks like 

( N-l) @xA z ) + {N-l)</i x ^(z)-NV'(z) + —^— + —^— + (regular terms of order 0(N°)) = 

z — A z — p 

(3.25) 

At this level of accuracy f\ ifl = f, the correction to / being a regular term of order 1/N. 
Then the correction — a = G\,n — G is the same for a or G, and the analysis of eq. (3.25) 
becomes perfectly identical to what was done in section 2 with eq. (2.17). We immediately 
write the differential equations that we obtain: 



1-logZ^ = eNinp(X)+NG(X) - ( 1 - {J^EMEd 

— log 2^ = e'Ninp(p) + N fa) - - 1 - fj, J= 



(3.26) 



Again it is clear that only the saddle-points (±, =f) need to be considered because the saddle- 
points (±, ±) are suppressed by a factor of 1/N . It is now convenient to introduce a modified 
kernel; noting that transformations: 

K(X,p)= 1 (X)K(X,p)-^- (3.27) 

do not affect values of determinants of type (1.2), we choose here 7(A) = exp^V^A)) so that 

K(X,p) ~ e- NV{fl) Z(X,p). (3.28) 
Using the saddle point equation (3.21), we can rewrite the differential equations: 



9 ■ log K± tT = ±Ninp(X) + NG(X)- ' 



(3.29) 



dX rv 7 ^ v 7 X-p 

^- log ^±, T = TNinp(p) - NG(X) - -i-^. 

Imposing the condition that ^(A, A) = Np(X), we finally get: 

R = sin(Nn(X-p)p(X)) em ^ mxl 

7T(A - p) 

This formula is a generalization of what was obtained in section 2 (which is the case A = 0). 
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The new factor (A — /i) j£r(A) can be absorbed in a redefinition of the kernel of the type (3.27). 



3.5. Multi-cut generalization. 

So far we have always assumed that the resolvent G has a single cut, i.e. the saddle point 
density of eigenvalues should be non-zero on a single interval Intuitively it seems clear 

that removing the "single cut" hypothesis, which is a long distance effect, should not change 
the short distance behavior of the kernel. Indeed, it has been checked in [19] that the short 
distance universality is preserved in the quartic {7(iV)-invariant case with two symmetric cuts 
(but the long distance behavior is modified). Here we shall argue that, in the more general 
models we consider, this change is irrelevant for short distance asymptotics. 

Let us start from eq. (3.25), which was derived without any assumption on the cuts of G 
or a. We shall study the correction C^(z) which satisfies 

W) = ^ (3.3D 

and its effect on the differential equations (the same line of reasoning will apply to C\{z)). 
Since our interest lies in the region A — ji ~ l/N, we can assume that both A and ji are in 
the same interval which is but one of the cuts of G. Of course A (resp. /i) must be 

slightly shifted in the imaginary direction to remove ambiguities, with a shift of ei (resp. e'i). 
As already noted, C^{z) cannot have a pole at z = n on the physical sheet. However, 

G^z) = \{C(z + iO) + C\z + iO)) = \{C{z - iO) + C*(z - iO)) (3.32) 

where C*(z) is the function connected to C(z) by the cut [a,(3] (note that the definition of 
C* does not depend on the cut chosen, because of the saddle point equation); so (3.31) can 
be continued to complex values of z and one concludes that C*(z) has a pole at z — [/,, with a 
residue of 1. 

It is now easy to derive the simplified differential equations for K e e > in an even more 
schematic way than before. If e = e', A will never get close to the pole of C* (figure 2); 
therefore there will be no 1/(A — //) term in the different equations and K e e > will of order 
0(N°), i.e. negligible. On the other hand, if e = — e', A will reach the pole of C* (figure 3); 
this time we obtain equations (3.29), and the usual short distance universality. 
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Fig. 2: Analytic structure of and a\ :fi when A and ji are on the same side of 
the cut When A ~ /i and the A on the physical sheet (indicated by a small 

arrow) approaches the cut, it does not reach the pole at z — /i on the other sheet. 




Fig. 3: Analytic structure of G\ :fi and ^ when A and /i are on opposite sides of 
the cut. This time the A on the physical sheet, as it approaches the real axis, "sees" 
the pole at z = 11 on the other sheet (through the cut). 

4. Multi-matrix generalization. 



The results obtained in the previous sections can be generalized to models of a chain of 
matrices. We shall only briefly discuss these, since they are less interesting physically. 
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4.1. The U(N) -invariant case. 

The first model considered is defined by the measure (with arbitrary potentials V^): 

M / / M M-l \ \ 

Yl d N2 M^ exp -Ntr - V^ m \M^) + M^M (m+ V . (4.1) 



It has a global U(N) invariance, and just as in section 1, can be treated by introducing the 
appropriate orthogonal polynomials. Here they are biorthogonal polynomials and Qk with 
respect to a non-local measure: 

. M 

I J] rfA(-)p fc (A( 1 ))gKA (M) )e^^^ (m) ( A(m) ) + ^^- A(ra)A ™ = 6 kl . (4.2) 

J m=l 

We can now define the distributions p n (Ai, ■ ■ ■ , A n ) of the eigenvalues of It is important to 
note that these distributions are different from the distributions of section 3, after (quenched) 

averaging over the external field A. They satisfy the usual determinant formulae (1.2) with 
the kernel 

/JV-i \ / M N N 

•* \i=l J \m=2i=l / 

A(X?\^ x det(eMN\P\f%^ (4.3) 

'M-\ 

fmK(m+l)^ \ An (M), 



H det(exp(iVAr ) Af +1) )) A(A^). 



We have expressed K directly as an integral over eigenvalues. Note that at this stage we can 
replace the determinants det(exp(iV Aj m) Aj m+1) )) in (4.3) with exp(iV £(Aj m) Aj m+1) )), which 
would lead to an expression of K in terms of biorthogonal polynomials: 



N-l 



K(\, fj.) = e-W»)) PkWPM (4.4) 



where 

M 



PM = I II ^(-)g fc (AM)e^ (2) ^^^ (m)A(m+1) ^^^ (ra) ) (4.5) 

m=2 

Here, as in section 3, we do not follow this path: we keep the determinants, in order to have 
saddle points. We then write down differential equations for the kernel, which can be simplified 
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in the region A — n ~ 1/N: the resulting equations are identical to eq. (3.29), and we find the 
same asymptotic expression (3.30). 

One could also consider correlations of eigenvalues of a matrix Af( m ) somewhere inside 
the chain, with presumably the same techniques and results. However, this has not been 
investigated in detail yet. 

4.2. The chain of matrices with an external field at one end. 

In this model we add to the measure (4.1) of the preceding section an external source term 
for MW: 

M / / M M-l \ \ 

Jl d N2 M^ exp -JVtr I - ^ V {m) (M^) + AM® + ^ M^M^ m+1 ^ . (4.6) 

m=l V \ m=l m—1 / J 

This measure is no more [/(A^)-invariant. We define again the distributions of eigenvalues of 
M^ l \ which satisfy determinant formulae (1.2): we give without proof the kernel 

„ /N-l \ / M N 

J \i=l J \m=2i=l 

det(exp(ATaAf ) )) A a) =A det(exp(iVA l (1) Af ) )) A (i)^ (4.7) 

M-l \ 

J] det(exp(iVA i (m) Af +1) )) A(A 4 (M) ). 

The differential equations will look a little more general; for example the equivalent of eq. 
(3.26) will be: 

^■]ogZ(X^) = (N-l)a^(X) 
— \w7A\. n.) = (N -1)> (2) 



-logZ(A,/i) = (iV-l)A A >). 



where 2 functions a(A*- 1 -*) and A*- 2 -*(A'- 1 -') must be introduced (the logarithmic derivatives of 
det (exp (iVaj A j 1 '')) and det(exp(iVA^ 1 ' ) Aj 2 ' ) ))). The rest of the analysis is the same. At the 
end, one redefines the kernel K with a transformation of type (3.27), using the saddle point 
equation: 

4(z) + ^ 2 \z) = V'(z) 
so that it satisfies the universal property (2.28). 
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5. Conclusion. 



Let us summarize our results. The model considered is that of a matrix coupled to an 
external field A; the latter has a smooth large N limit characterized by a limiting density of 
eigenvalues. A kernel K(X,/j,) is defined such that the determinant formulae (1.2) hold. In the 
simplest case where A — 0, the large N form of K(X,fi) for all A and /x, first found in [11], is 
reproduced here. In the case of a non-zero A, we restrict ourselves to the region A — /i ~ 1/N 
(it is not clear that, for general A and V, the long distance behavior of K, even after smoothing 
the oscillations, should be interesting, e.g. should exhibit any kind of universal behavior). The 
asymptotic form (3.30) is obtained, extending the level spacing universality to this class of 
models. 

The key ingredient of the derivation of the short distance universality is of course the 
existence of the kernel. But it seems reasonable to assume that the level-spacing universality 
(observed experimentally for a broad range of systems) should be true for very general matrix 
models, in which the correlation functions do not satisfy (1.2). The problem is then to manage 
to compute the level-spacing distribution P(s), even though it is no more simply related 
(through K, cf eq. (1.5)) to the correlation functions, that is the naturally calculable quantities 
of the model. This suggests that a totally different approach is probably necessary. Here, let 
us mention that the question of knowing how far the universality of P(s) can be extended 
is reminiscent of the question of knowing what is the domain of attraction of a fixed point 
in renormalization group theory. What we have shown is that in our matrix models, both 
non-gaussian terms and terms explicitly breaking the ?7(iV)-invariance are irrelevant in the 
large N limit and lead to the gaussian £/(iV)-invariant fixed point. Maybe RG methods can 
be applied here too (cf [9]) and allow a much more general approach to the problem. 

In the case of multi-matrix models, the same short distance behavior is found for the 
correlation functions of the eigenvalues of a given matrix, here the first matrix in a chain of 
matrices. Of course, it would be interesting to investigate the more general problem of the 
correlations between eigenvalues of different matrices in the chain. This has been already done 
in the gaussian case [17]. The conclusion is that the interesting limit is that of an infinite 
chain, which tends to the c = 1 matrix model 



Z 




21 



One should then study the correlations of the eigenvalues of M(t\) and M(t2) as t\ and t<i are 
close to each other. In order to apply the methods used in this paper, one will of course have 
to generalize the determinant formulae (1.2) to include eigenvalues of different matrices. 
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APPENDIX 1. Analytic structure and functional inverses. 

Let us define the function a(A) used in section 3, and its functional inverse A(a). We shall 
assume, as in [13], that the logarithm of the Itzykson-Zuber integral 

W = —\og det{eMNaiX * )) (All) 
AT2 10g A(A,)A(a z ) ^> 

has a smooth large N limit, so that it depends only on the distributions of eigenvalues p(A) 
and p(a) of M and A. Then one can define functional derivatives of W with respect to p and 
p. It is clear (eq. (3.16)) that 

« A > " m = Txlk) WM (AL2) 

The r. h. s. can be extended to complex values of A and has no cut on [a, f3] (the support of 
the density p), therefore we can define a(A) by removing the slashes in (Al.2). Note that this 
corresponds to the very simple finite N definition: 

a(A) = ^^logdet(exp(A^A;))Aiv + i=A (A1.3) 

The only problem of this definition is that there are now N eigenvalues a\ and N+ 1 eigenvalues 
\f one should add one a\ or remove one Aj. In the large N limit, this problems disappears, 
since if we believe that W depends only on p(A) and p(a), then one can redefine the eigenvalues 
a\ (or Aj) to add an eigenvalue a\ (or remove one Aj), keeping the densities p(a) and p(A) fixed. 
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All that has been done so far is symmetric in the exchange of A and M, so we can also 
define A (a): 

X ^=^ + Talk) WM <AL4) 

where G(a) is the resolvent of A, and the functional derivative has been again extended to 
complex values of a. 

Let us now discuss the analytical structure of a(A) and A(a). On [a, (3], according to (Al.2), 
a(A) has the same cut as G(A), i.e. 

Ima(A±iO) = ±7rp(A). (A1.5) 

Likewise, A(a) has the same cut as G(a). We can then define a*(A) (resp. A*(a)), to be the 
functions on the other side of the cut of G (resp . G). If we now consider the matrix model 
with an external source term (measure (3.1)), a* (A) satisfies an additional constraint which is 
the saddle point equation (3.22). One can see that this implies that a*(A), just like G(A), has 
no other cut in the whole complex plane than that on [a, (3}. On the other hand, A* (a) is not 
constrained and therefore it may have more cuts (leading to other sheets). 

Finally, it can be shown by studying the large N limit of the Itzykson-Zuber integral 
that A(a) and a(A) (as multi-valued functions) are functional inverses of each other. This 
can be thought of as a generalization of the inversion relation found in [13], even though the 
connection is non-trivial. Here we shall derive this relation in an elementary fashion. 

We rewrite definition (Al.3) explicitly: 

a(A) = E ^ + r ( A1 -6) 

where £^+1 is the determinant det(exp(Aa^Aj)) with 1 < k < N + 1, k ^ I and 1 < i < N. 
Now eq. (Al.l) can be applied to E^n+i: 

>^^M^r WM+0[llN2) (AL7) 
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where p\ is the density of the k ^ I; computing the 1/N correction to log-E^jv+l one finds 
\ogE ljN+1 = C- N-^—W[p,~p] ~ V log(a fc - a,) (A1.8) 

where C is independent of /. Finally, using the standard trick which is to replace the sum over 
/ with a contour integral, (Al.6) becomes 

SdaG(a)ae NaX e~ N rHa')da' 

a A = ^ — 77 —- A1.9 

1 ^ $ daG(a) e"<*e- N f X ^ da ' 

The saddle point equation gives simply: a(A) = a with A(a) = A, which proves the functional 
inversion relation. 

A last remark: the analytic structure described here can be made more explicit in the 
gaussian case (with an external field), by using Pastur's self-consistent relation [5]. 



APPENDIX 2. Connection with characters. 

One may notice the strong resemblance between equations discussed in section 3.3 and 
appendix 1 and large N character relations found by Kazakov et al. in [16]. To establish the 
connection, we shall now rederive in a very simple manner the latter relations. 

A representation of U(N) can be described by its highest weights mi > rri2 > . . . > m^. 
The corresponding character is defined by: 

(A-) 

where the Z{ are eigenvalues of the N x N matrix M , and the hi are the shifted weights 
hi = mi + N — i. 

If the weights are positive (mjv > 0, which we shall simply write h > 0) then we can use 
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the identity 

E Xh'(M) Xh '(Mo) = exp I f; - trM'tr M* ) (A2.2) 

h'>0 \q=l q 

and the orthogonality of characters to compute 



Xh(Mo) = J dMxh(M) I Xhi{M) Xh ,{M { 

MeU(N) \h'>0 

= J dMxh(M)exp(NtrV(M)) 

MeU(N) 



(A2.3) 



The potential V(z) = Y^ =l t q z q /q is related to M by t q = ^tr(M^). 

We now diagonalize M and integrate over the unitary group; we are left with the eigenvalues 
Zi of M, \zi\ = 1; using the definition (A2.1) we obtain: 

r N-l , 

XH(Mo) = I II det(*r™')e"£« V{Zi) . (A2.4) 

(the shifted weights /ij have been rescaled by a factor of AT). This is to be compared with 
partition function (3.15). The study of the two cases (and in particular the analytic structure 
that arises) being very similar, we shall now skip the details of the derivation. 

One first writes down a saddle point equation: 

$(z) + ${z) + V'{z) = (A2.5) 
where the function k(z) is defined by: 

^) = lAlogdet(^). (A2.6) 

k(z) has the same cut as G(z), i.e. where the saddle point density of eigenvalues is non-zero. 
We then have 

k*(z) + G(z) = V'(z) (A2.7) 
on the whole complex plane, (k* is the function on the other side of the cut of k). We make 
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the trivial redefinition h(z) = 



—zk*(z), and expand in powers of z: 



oo oo / 1 N—l 

9=1 g=1 V 1=0 ' (A2.8) 

OO OO ^ V ' 

= e + 1 + E z ~ 9 j^ q dT log Xft(Mo) 

which is precisely the result obtained by much more complicated methods in [16] (in their 
notation z = G -1 .) 

Finally, the same functional inversion argument as in appendix 1 applies here, proving that 
z(h) defined by z(h) = expw(h) with 

Mhi) = -±-£-]og X h(Mo) (A2.9) 
(cf [16] for a similar definition of G(h)) is the functional inverse of h(z). 
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